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I. INTRODUCTION 


The existence of cosmic rays with energies above the Greisen-Zatsepin-Kuzmin (GZK) cutoff [Ei is one of the 
most important open problems in high-energy astrophysics [ 1 , 1 ^. These ultra-high energy cosmic rays (UHECRs) 
may be a tantalizing hint of novel and very powerful astrophysical accelerators, or they may be harbingers of new 
microphysics via the decays of metastable supermassive particles. It is remarkable that we still do not know whether 
the UHECRs originate from macrophysics or microphysics. If there is no GZK cutoff, as suggested by the AGASA 
data Q, the sources of the UHECRs would need to be local. In this case, since local magnetic fields are unlikely 
to have deflected significantly their directions of propagation, the UHECRs would ‘remember’ the directions of their 
sources. Thus, one would expect some anisotropy in the arrival directions of the UHECRs, associated either with 
discrete energetic astrophysical sources nearby, such as BL Lac objects Q, or the distribution of (mainly galactic) 
superheavy dark matter. No significant anisotropy of the UHECRs has yet been seen, but the existing experiments 
have insufficient statistics to exclude one at the expected level On the other hand, the GZK cutoff may be present 
in the HiRes data [E, in which case no exotic microphysics may be needed, and any astrophysical sources would be 
less restricted and more difficult to trace. The first batch of Auger Q data are inconclusive on the possible existence 
of the GZK cutoff. 

Superheavy particles of the type required could have been produced gravitationally around the end of inflation . 
Particularly interesting candidates for such superheavy particles are ‘cryptons’, bound states of the fractionally- 
charged constituents that arise generically [l^ in the hidden sectors of models of particle physics derived from the 
heterotic string Cryptons arising in the hidden sector of a heterotic string-derived flipped SU(5) model may have 
exactly the right properties to play this role [ll|, [l^. These ‘flipped cryptons’ are bound by SU(4) hidden-sector 
interactions, and include four-constituent meta-stable bound states called tetrons that are analogues of the three- 
constituent baryons of QCD, as well as two-constituent meson-like states. Indeed, it was within this flipped SU(5) 
model that the confinement solution to avoiding the stringent experimental limits placed on fractional ch arg es was 
first pointed out, and this model remains the only example to have been worked through in any detail [Til. Il2j|. 

In general, tetrons may decay through A^‘^-order non-renormalizable operators in the superpotential, which would 
yield lifetimes that are expected to be of the order of 


2-N 

string 


{-) 

mx \mxJ 


2(N-3) 


( 1 ) 


where mx is the tetron mass and Ms ~ 10^® GeV is the string scale. The a-dependent factor reflects the expected 
dependence of high-order superpotential terms on the effective gauge coupling g. The mass scale associated with these 


Such states may also be a generic feature of models constructed from intersecting D-branes. 


1 




2 


states is estimated using the renormalization group for the SU(4) interactions to be A ^ 10^^ — 10^^ GeV, just in the 
right range for their decays to produce the UHECRs. 

The lifetimes of neutral tetrons without electric charge has been estimated to lie in the range tq ^ 10^^ —10^^ years, 
so that they may still be present in the Universe today, and might produce the necessary flux of UHECRs if they are 
sufficiently abundant. We have shown in our previous work [l^ that the mesons and charged tetrons - whose present- 
day abundances are subject to very stringent experimental limits - would have decayed with short lifetimes early in 
the Universe. In the course of studying the possible tetron lifetimes, we identified various 10*^-order superpotential 
operators that might govern neutral tetron decays. Thus, in this specific model of flipped cryptons we are able to 
go beyond generic statements regarding the injected UHECR particle spectra that may result from their decays, and 
make a number of specific predictions. 

This enables us to address an important experimental constraint on such crypton models of the UHECRs. Although 
‘top-down’ models such as crypton decays may appear to be natural explanations for the UHECRs (if they exist), they 
generically share a potential drawback. The spectra of UHECRs that they produce might be expected to have large 
photon fractions, in possible conflict with the observation that most of the UHECR primaries appear to be protons 
or nuclei. The Auger collaboration has recently set an upper limit of 26% on the photon fraction above 10^® eV [H, 
and the limit may even be as low as 7-14% [ij]), while an upper limit of 50% at energies above 10^® GeV has been 
set by the AGASA collaboration [l^. 

In this paper, we first give a review of relevant aspects of the flipped SU(5) heterotic string model. We then analyze 
the specific primary multi-body decay modes governed by the various different 10*^-order superpotential operators 
found in our previous paper. We then calculate the UHECR particle spectra that would be injected by these decays 
using one of the most detailed and complete codes currently available |16l | , paying particular attention to the photon 
fraction. The total UHECR spectra obtained from the various superpotential terms do not differ greatly. On the 
other hand, we find that different decay operators may give rather different photon fractions, particularly at the 
highest energies. However, in every case, we find that the calculated spectra after supersymmetric evolution and 
fragmentation are compatible with the published upper limits on the photon fractions in various energy ranges, when 
we include the UHECR background resulting from a homogenous extragalactic distribution of sources and incorporate 
the pile-up expected from the GZK effect. There is no need to appeal to the cosmic radio background to absorb a 
significant fraction of the photons in order to bring them below the AGASA and other limits. 


II. GENERIC SUPER-HEAVY RELIC DECAY 

The basic idea in generic ‘top-down’ explanations (see [EEl) is that the UHECR are produced via the decay of 
some relic particles or topological defects left over from the inflationary epoch and which are locally clustered in the 
galactic halo as cold, dark matter with an over-density ~ 10'^“®. The lifetime of such relics must exceed the 

present age of the universe in order for them to exist today in sufficient abundance, however the lifetime must not be 
too large so that the decay rates produced are too small to produce the UHECR. Furthermore, the relic mass must 
be at least Mx > 10^® GeV in order to produce the UHECR energies observed. Typically, the lifetime of a particle 
is expected to be inversely proportional to it’s mass, r ^ 1/M. Clearly it is not easy to have a particle with both a 
large enough mass and a decay lifetime in the right range to produce the UHECR. However, as pointed out in the 
Introduction, flipped SU{5) cryptons satisfy both of these criteria, which makes them very attractive as a top-down 
explanation of the UHECR. Indeed the ‘flipped’ crypton is probably the most natural and most physically motivated 
of any top-down candidate. 

There are three generic statements that can be presently made about a decaying super-heavy X particle explanation 
for the UHECR: 

1. Since the super-heavy relics may accumulate locally in the galactic halo with an over-density 10^ ® over the 
cosmological average, they may avoid the GZK cutoff. 

2. Due to the displacement of our solar system with respect to the galactic plane, there should be some anisotropy 
in the arrival directions of the UHECR with respect to the galactic center. 

3. Photons tend to be the dominant component of the UHECR flux produced by the super-heavy X decay. However, 
the photons may scatter off the galactic radio background, which is poorly measured, and thus may be somewhat 
attenuated. 

The injection spectrum produced by such a decaying super-heavy relic V-particles with number density nx and 
lifetime tx is proportional to the inclusive decay width: 

nx 1 dr{X ^ gi-\ -) 
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If a spherical halo of radius Rhaio and uniform number density nx is assumed, then the galactic halo contribution to 
the UHECR will be given by 

jhalo ^ ( 3 ) 

47r 

In general, the X-particles will decay into one or more partons which hadronize into other particles g of the MSSM, 
which carry a fraction x of the maximum available momentum Mx /2 and a fraction z of the parton momentum. For 
such a decay, the inclusive decay width can be factored as 
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where is the fragmentation function (FF) for particles of type g into particles of type a, and g, is the energy 

scale, most appropriately taken to be equal to the X particle mass, g = Mx- The evolution of the fragmentation 
function is governed by the DGLAP equations, which may be extended to include the MSSM. Thus, the determination 
of the expected UHFCR flux from the super-heavy X decay essentially becomes the problem of starting with a set of 
initial decay partons and evolving the decay cascades via the fragmentation functions to find the end decay products 
and energy distribution. To evolve the fragmentation functions up to the energy of the super-heavy X decay, the 
DGLAP equations must be solved numerically. Several grou ps have done such calculations for generic initial decay 
partons (usually into a quark-antiquark pair) [^■ l22j|. Perhaps the best such code is SHdecay [I^. This code 

calculates the fragmentation into the seven stable MSSM particles (p, 7 , e, neutralino LSP for any given 

initial decay parton. In the case of flipped SU{5) cryptons, we have a specific model where the initial decay partons 
in the cascade are known. 

In addition to the UHFCR flux produced by the super-heavy decay from X particles clustered in our galactic halo, 
there may be a background flux from sources outside of our galaxy, perhaps super-heavy X decay in other galactic 
halos or in intergalactic regions. Generally, this flux is assumed to be due to a homogenous distribution of sources 
and exhibits a characteristic GZK pileup due to the fact that they are produced non-locally [s^. Since the GZK 
attenuation is much more severe for photons than for nucleons, this background should be comprised primarily of 
nucleons. Thus, in this scenario the UHFCR flux observed on Farth will be the sum of this extragalactic background 
and the local galactic flux from decaying relics clustered in our galactic halo. Due to the extragalactic component, 
super-heavy X particle decay may only unambiguously explain the UHECR flux for energies E > 4-10^® eV. However, 
we note that the extragalactic component may also be due partially to non-local super-heavy X decay, as well from 
astrophysical sources (‘bottom-up’ production). Thus, it is more accurate to say that a distinct signal of super-heavy 
X decay within our galactic halo would be the existence of an excess of events above this energy. The lack of any 
events above this energy would not rule out the presence of a top-down component, but it would not provide an 
unambiguous reason for the introduction of such a mechanism. 


III. STRING-DERIVED FLIPPED SU(5) 


Even before the first string models, there was strong interest in flipped SU(5) as a possible grand unified theory, 
primarily because it did not require large Higgs representations and avoided the constraints of minimal SU(5) (23| 
without the extra gauge interactions required in larger groups such as SO(IO) [H, [1^. Interest in flipped SU(5) in¬ 
creased within the context of string theory, since simple string constructions could not provide the adjoint and larger 
Higgs representations required by other grand unified theories, whereas they could provide the Higgs representations 
required for flipped SU(5). It was, moreover, observed that flipped SU(5) provided a natural ‘missing-partner’ mech¬ 
anism for splitting the electroweak-doublet and colour-triplet fields in its five-dimensional Higgs representations [ 2 ^. 
Heterotic string-derived flipped SU(5) remains among the most fully developed and realistic models yet derived from 
string. 

Flipped SU(5) derives its name from the flipping of the quark and lepton assignments relative to those in the 
minimal SU(5) GUT [ 2 ^: ^ dL,d‘i and eL,e£ for all three generations. This requires the gauge 

group to become SU(5) x U(l), so as to accomodate the charges of the quarks and leptons, which fill up the 5, 10 
and 1 representations of SU(5): 


/s — 



e 

Me) L 



II = e\ 


L’ 


( 5 ) 







4 


The presence of a neutral component in the 10 allows spontaneous GUT symmetry breaking to be achieved using a 10 
and a 10 of superheavy Higgs fields with the same and opposite U(l) hypercharges. Their neutral components develop 
large vacuum expectation values (few’s) breaking SU(5) x U(l) ^ SU(3) x SU(2) x U(l), while the electroweak 
spontaneous breaking occurs through the Higgs doublets H 2 and H^: 

h5 = {H2,H3} ; h5 = {H2,H3}. (6) 

The resulting economical doublet-triplet splitting mechanism gives a large mass to the Higgs triplets (HsjH^) by 
coupling them to states in 10 and 10 Higgs representations, through trilinear superpotential couplings of the form 

FFh ^ d’H H 3 (7) 

FFh^dj,{D%)H-^, (8) 

while keeping light the Higgs doublets (H 2 ,H 2 ) light. 

The absence of any mixing between the Higgs triplets (H 3 , H 3 ) in this dynamic doublet-triplet splitting mechanism 
provides a natural suppression of the d = 5 operators that might otherwise mediate rapid proton decay, so that 
flipped SU(5) is probably the simplest unified model that can satisfy the experimental limits placed on the proton 
lifetime |^ . 

String-derived flipped SU(5) belongs to a class of models constructed using free fermions on the world sheet, 
corres pon ding to compactification on the Z 2 x Z 2 orbifold at the maximally-symmetric point in the Narain moduli 
space |28l |. Although this model was originally constructed in the weak-coupling limit, it is possible that it may be 
elevated in the strong-coupling limit to an authentic M-theory model and may at some point make contact with 
models based on D-brane constructions. 

The full gauge group of the model is SU(5) x U(l) x U(l)'*x SO(IO) x SO( 6 ) SU(4)], where the latter two 

factors are confining hidden-sector gauge interactions. The matter spectrum comprises the following fields: 

(i) Observable sector: The conventional matter fields may be regarded as three 16 representations of SO(IO) that 
D SU(5) X U(l) chiral multiplets Ui(10, i), /j(5,— |), ^(’(1, |)(i = 1,2,3); extra matter fields ^ 4 ( 10 ,^), / 4 ( 5 , |), 
[4(1, — |) and ^5(10, —5), 75(5, —|), ^5(1, |); and four Higgs-like fields in the 10 representation of SO(IO), that D 
/ii(5, —1), /ii(5, 1), z = 1, 2, 3,45. In our realization of the model, we make the following flavour identifications of the 
Standard Model states with the various string representations: 
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where /3 indicates a mixture of fields with the indices 1 and 3. The light Higgs doublets hd and /i„ (where hd is the 
Higgs doublet responsible for giving mass to the down-type quarks and leptons and gives mass to the up-type 
quarks) are contained in in the hi and (145 string pentaplet representations, respectively. 

(ii) Singlets: There are ten gauge-singlet fields (/) 45 , (/>+, (j)i{i = 1,2,3,4), $ 12 , $ 23 , $ 31 , their ten ‘barred’ 
counterparts, and five extra fields $/(/ = 1 ■ • • 5). 

(iii) Hidden sector: This contains 22 matter fields in the following representations of S'O(10)zi ® SU{4)h'- 7i(10,l), 
Ai(l, 6)(z = 1 • • • 5); Fi{l, 4), Fi{l, 4){i = 1 ■ • • 6). Flat potential directions along which the anomalous combination 
of hypercharges U(1 )a is cancelled induce masses that are generally near the string scale for some, but not all, of 
these states. Depending upon the number of Ti and states remaining massless, the SO(IO) condensate scale is 
2^014-15 QgY 2he SU(4) condensate scale is 10^^“^^ GeV [ 2 ^. The 73^5 and Fs^^ states always remain massless 
down to the condensate scale. The U(l)i charges and hypercharge assignments are shown in the Table below. 
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Table: T/ie spectrum of hidden matter fields that are massless at the string scale in the revamped flipped SU(5) 
model. We display the quantum numbers under the hidden gauge group SO(IO) x S0(6) x U(lfl, and subscripts 
indicate the electric charges. 


The hidden-sctor matter fields are confined into crypton bound states. These occur in ‘cryptospin’ multiplets with 
different permutations of confined constituents, analogous to the flavour SU(2), SU(3) and SU(4) multiplets of bound 


states in QCD. The cryptospin multiplets contain doubly-charged tetrons 



4' — F 3 F 3 P 3 P 3 , 

4-++ = F 5 F 5 F 3 F 5 , 

(10) 


^++ = F 3 F 3 F 3 F 3 , 

4" = hhhh, 

(11) 

and singly-charged tetrons 

4'^ = F 3 F 3 F 3 F 3 , 

4' = F 3 F 3 F 3 F 3 , 

(12) 


4 = F 3 F 3 F 3 F 3 , 

4'’" = F 3 F 3 F 3 F 3 , 

(13) 

as well as neutral tetrons 

^0 = F 3 F 3 F 3 F 3 , 

40 = p3p3hh- 

(14) 

We have shown previously that the charged tetrons could have decayed with rather short lifetimes early in the history 
of the Universe |l^. and so avoid any problems with charged dark matter particles [s^. 

However, the neutral tetrons can decay only via higher-order operators in the superpotential, which we have 
identified in [T^, and may have lifetimes exceeding the age of the Universe. This may make them good candidates 
for cold dark matter, and their slow decays are a possible source of the UHECRs. 

In the next Section, we turn our attention to these interactions and examine the UHECR energy spectra that may 


result from these specific tetron decays. 


IV. UHECR INJECTION SPECTRA 

We have previously found [l2| the following 10*^-order superpotential operators through which the neutral tetrons 
may decay: 

4'°[F2F2$3i045(^ hi + F 2 F 2 ^ 2 ‘ifiibfi^ hi + T2^3-p3'/’4 045/2 + 

^ 44 ’ 230450 ~/ 45/5 + (^ 310450 ” 

+ 4>23 0450''’)^i(/2^2 + h^V) + ‘J’23 0450 


( 15 ) 
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/li + F2i^2^23'^45</^'^^1 + ^ 2^20 / 11 / 11/145 + (16) 

F4F4$3i(/)45())“'"/ii + FiFi(j)^hihihi^ + -F4‘53i(()45(()^^45/5 + 

Fi(j)^hihi^h/^^J^ + F^^ hihihil^ + 

{^3l(j)45'P~^hi + (I)^hihihi^)fil1 + ($ 31 ()) 45 (^ hi + $ 230450~''^1 + 

(/)■'■/ 145 / 145/145 + 0 /ll/ll/l45)(/2^2 + /s^s)]- 

Using the flavour identifications we outlined above, these operators would give rise to the following neutral tetron 
decay inodes: 

^ T hd 5'° ^ e/^ hd 5'° ^ hd (f)^, (17) 

6 /id hd hu4>, 5'° ^ hu ^tFhu h^ hd 4>, 

5'° ^ c c'^ d III® —> s hd 4)^. 


We note that there are several different possible decay modes, any of which may be dominant, depending on unknown 
features of the model dynamics that determine the relative values of their coefficients. In particular, the most 
important tetron decays could be into either leptons or quarks, and there are many different possibilities for the 
dominant flavours. 

We plot in Figs. 1 to 8 below the expected UHECR energy spectra of photons and nucleons due to each of these 
possible tetron decay modes, as well as the maximum photon fractions expected. The energy spectra were calculated 
for a mass Mx = 2 • 10^^ GeV, using the fragmentation functions D‘^{x,M]^) generated by the code SHdecay 0 - 
This code calculates the fragmentation into the seven stable MSSM particles {p, 7 , e, neutralino LSP y, Ve, Vt) 
for any given initial decay parton. The many-body decays distribute the total decay energy Mx among the different 
particles. We include Higgs decays, but we ignore the decays of the singlet fields, except to take into account their 
kinematical effects on the primary quark and lepton spectra. We follow in estimating the probability density 
Pniz) that one decay parton carries off a fraction z of the total available decay energy Mx- 

p„(z) = (n - l)(n - 2)z(l - z)”“^ (18) 


for n > 3 decay partons. The resulting flux from the emission of a given decay parton is then 

E^J\E)=Bx^ (19) 

where i = {jp, 7 , e, y, Ve, Vr)- 

To obtain the total UHECR spectrum, we add to this the background flux of nucleons that would be expected to 
result from a homogenous distribution of extragalactic sources that exhibits the distinctive pile-up due to the GZK 
effect [s^ The constant B in (I19p is a normalization coefficent determined by the tetron number density and 
lifetime, viz 

( 20 ) 


This dimensional coefficient B is not determined a priori, and must be fitted to the experimental data. In each of 
Figs. 1 to 8 , we show the total spectrum obtained by summing the background and the fluxes of nucleons and photons 
resulting from tetron decay, and in a second panel we display the gamma fractions: 7/(7 +p). We have assumed no 
photon attentuation in the calculated spectra, although a strong attenuation cannot be excluded IM, because the 
galactic radio background has never been accurately measured and its intensity is largely unknown [35j |. 

Figs. 1 to 6 are for quark primaries, and are ordered according to the masses of the quarks involved. In the case of 
the b quark, we show in Figs. 3 and 4 plots for two superpotential operators with different numbers of accompanying 
Higgs fields: the two plots are rather similar, and the same is true of the two plots shown in Figs. 5 and 6 for primary 
t quarks. We are thus led to hope that including the (model-dependent) decays of the singlets <p would not have large 
effects. The plots for lepton primaries shown in Figs. 7 and 8 are more distinctive, in that the photon fractions rise 
to much larger values at energies above 10^° eV 


^ However, we note that this model is likely to come under pressure from upper limits on high-energy cosmic-ray neutrinos |,S3ll - private 
communication from S. Sarkar, see also http;//www-thphys.physics.ox.ac.uk/users/SubirSeirkar/talks/niunich05.pdf 
® The photon fractions for second-generation quark primaries look somewhat flatter than those for b and t quarks, but this difference is 
probably within the modelling uncertainties. 
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FIG. 1: The top panel shows the total UHECR spectrum 
and the bottom panel the photon fraction from the decay 
mode —> s hd ■ 


FIG. 2: The top panel shows the total UHEGR spectrum 
and the bottom panel the photon fraction from the decay 
mode ^ c d (F (j)^. 



FIG. 3: The top panel shows the total UHEGR spectrum 
and the bottom panel the photon fraction from the decay 
mode hd (j>^ ■ 



FIG. 4: The top panel shows the total UHEGR spectrum 
and the bottom panel the photon fraction from the decay 
mode iP® —> fe hd hd hu4>- 





























FIG. 5: The top panel shows the total UHECR spectrum 
and the bottom panel the photon fraction from the decay 
mode —> t 0®. 


FIG. 6: The top panel shows the total UHEGR spectrum 
and the bottom panel the photon fraction from the decay 
mode 'h® —> t hu hu hd 4>- 




FIG. 7: The top panel shows the total UHECR spectrum 
and the bottom panel the photon fraction from the decay 
mode 4'° —> e//r hd ■ 


FIG. 8: The top panel shows the total UHECR spectrum 
and the bottom panel the photon fraction from the decay 
mode ^ T hd </>^. 
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FIG. 9: A comparison with the available data on the UHECRs from the Fly’s 
Eye, HiRes, AGAS and Auger experiments with the crypton decay model —» 
b fo'’ ha for Mx = 10^® GeV. 


In Fig. 9 we compare the spectrum for one of the operators with primary b quarks, calculated for a crjmton mass 
of 10^^ GeV, with experimental data from the Fly’s Eye, HiRes, AGASA, and Auger experiments 011 Si- The 
AGASA flux has been scaled by a factor of 0.55 for consistency with the other data, and the normalizations for the 
crypton decay contributions to these spectra has been adjusted for the different crypton masses. The limited statistics 
for UHECRs with energies > 10^® eV available in the present data sets do not offer any clear discrimination between 
crypton masses in the range 2 x 10^^ GeV > Mx > 10^^ GeV. In the case of a crypton mass ~ 10^^ GeV there is 
no clear signal of a crypton contribution to the UHECR since the flux from such a decay is essentially buried within 
the background from homogenous extragalatic sources. A clear signal of crypton decay, at least in this model, would 
require a lower limit on the crypton mass Mx > 5 • 10^® GeV in order to provide an excess of events above 4 • 10^® eV 
that could not be attributable to extragalactic astrophysical acceleration mechanisms. 


V. CONCLUSIONS 

We have carried as far as is possible at present the modelling of flipped crypton decay contributions to UHECRs, 
including all the possible 10‘^-order superpotential operators. The experimental data presently available are consistent 
with all the decay modes possible in this crypton framework. The total UHECR spectra are consistent with a 
contribution from cryptons weighing between 2 x 10^^ GeV and 10^® GeV, although only a crypton mass Mx > 
5 • 10^® GeV would provide an unambiguous signal over conventional explanations. The available upper limits on the 
possible photon fraction do not exclude any of the crypton models we have studied. 

In the future, the larger data set expected from Auger may be able to discriminate between crypton decays and 
other models of UHECRs, and also among different crypton models themselves. Greater statistics will enable the 
UHECR anisotropy to be measured with sufficient accuracy to discriminate crypton decay from a uniform distribution 
of astrophysical sources, and more accurate measurements of the photon fraction at higher energies might offer some 
discrimination between models with lepton and quark primaries, as seen by comparing Figs. I to 6 with Figs. 7 and 
8 above. 

Thus there is hope that, in the near future, we may finally learn whether UHECRs have a macrophysical origin or 
a microphysical origin and, in the latter case, may start to discriminate between different microphysical models. 
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